Mixing by polymers: experimental test of decay regime of mixing. 
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By using high molecular weight fluorescent passive tracers with different diffusion coefficients and 
by changing the fluid velocity we study dependence of a characteristic mixing length on the Peclet 
number, Pe, which controls the mixing efliciency. The mixing length is found to be related to Pe 
by a power law, Lmix oc pg0.26±o.oi ^ increases faster than expected for an unbounded chaotic 
flow. Role of the boundaries in the mixing length abnormal growth is clarified. The experimental 
findings are in a good quantitative agreement with the recent theoretical predictions. 



PACS numbers: 47.27.-1,47.50. +d,83.50.-v 



Mixing of a low difFusivity passive tracer by a turbulent 
flow, particularly in a region of small scales, attracted 
recently significant attention from both theory [1-5] and 
experiment [6,7]. This so-called Batchelor regime of mix- 
ing exists in a viscous-convective range of wave num- 
bers, kx < k < ks, where the random-in-time veloc- 
ity field can be treated as spatially homogeneous. Here 
kK = (e/j/3)i/4 and kg = {e/uD^/^ are the Kol- 
mogorov dissipation and the Batchelor wave numbers, 
respectively, v and D are the kinematic viscosity and the 
diffusion constant, respectively, and e is the rate of en- 
ergy dissipation [8]. The inverse injection scale should 
be used here instead of fc^ , if the former is larger than 
kK- The range of the Batchelor regime is defined by 
the value of the Schmidt number Sc = I'/D, which can 
be rather large for a viscous fluid with a low diffusivity 
tracer. The smoothness of the velocity field for devel- 
oped turbulence below the Kolmogorov scale is explained 
by an exponential decay of the velocity spectrum there. 
The Batchelor-Kraichnan approach [8] can be extended 
to an algebraic type velocity spectrum cx fc~* with (5 > 3. 
Then fiuctuation spectrum of the velocity gradients scale 
as k^^^^, so that the flow becomes increasingly homo- 
geneous on small scales, and the mixing occurs mainly 
due to eddies corresponding to the border of the velocity 
smoothness. Batchelor and Kraichnan [8] considered a 
stationary regime of turbulent mixing with a tracer con- 
tinuously injected into a turbulent flow on a large scale. 
Then the mixing is produced by stretching and folding of 
a tracer blob by the random velocity field creating com- 
ponents of higher and higher wave numbers up to fc^. 
A more recent theory of the passive scalar decay in the 
Batchelor regime of mixing, which is more relevant to a 
real experiment, makes predictions about the statistics 
of the passive scalar [5] . 

Various theoretical predictions for the stationary and 
decay cases of the Batchelor regime of mixing were veri- 
fied recently in macro-systems [6,7]. However, the theo- 
ries mentioned [8,2,5] consider an unbounded case only. 



As shown just very recently [9], a non-uniform velocity 
distribution particularly close to a wall, can alter signifi- 
cantly the efficiency of mixing in the decay regime due to 
turbulent advection. Indeed, due to a reduced velocity 
near the wall the boundary layer becomes a sink for the 
passive tracer. This excess of the tracer is intermittently 
injected from the boundary layer into the bulk. Thus, 
the tracer decay, e.g. along a channel, is controlled by 
the rate of tracer injection from the boundary layer that 
causes a significant slow down of the decay compared 
with an unbounded case. As a result, dependence of 
the mixing length on the Peclet number, Pe = Vayd/ D, 
changes from logarithmic to algebraic, Lmix c>c Pe-^^'^. 
Here Pe defines the relative rate of advection versus dif- 
fusive transport in a flow, Vav is the average longitudinal 
velocity, and d is the channel width. Besides, the bound- 
ary layer width for mixing in the decay regime in scales 
as oc Pe"^/**. In order to verify these theoretical pre- 
dictions for scaling of the mixing length, of the width of 
the boundary layer, and of the correlation function with 
Pe, we used a regime of elastic turbulence [10], which 
occurs at very low Reynolds numbers. Re ^ 1, in a vis- 
cous fluid containing a small addition of high molecular 
weight polymers, similar to the experiment reported in 
Ref. [7]. 

Solutions of flexible polymers are viscoelastic fluids [11]. 
The distinctive characteristic of such fluids is the pres- 
ence of elastic stresses in a flow that relax with some 
macroscopic characteristic time and increase non-linearly 
with the flow rate [11]. The latter is the result of the 
nonlinear mechanical relation between the elastic stresses 
and the rate of deformation of polymer molecules in the 
flow. As shown theoretically and experimentally, the 
non-linear elasticity of the polymer solution can lead 
to a purely elastic flow instability [12,13], and even to 
elastically-driven turbulence as discovered recently [10]. 
These nonlinear elastic effects are observable at an ar- 
bitrarily small Re, being practically independent of it, 
but strongly depend on Weissenberg number, Wi = A- 7. 
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The latter reflects the competition between non-linear 
elastic stresses and viscous stresses for the polymer part. 

Here A is the polymer relaxation time, and 7 is the shear 
rate. It was also shown recently [7] that a small addition 
of polymer to a viscous fluid can lead to a very eflicient 
mixing in a regime of the elastic turbulence generated in 
a curvilinear macro-channel, and that this system hap- 
pens to be an ideal test case for the Batchelor regime of 
mixing. 

In order to drastically reduce time for preparation of 
polymer solutions and at the same time to reduce the 
amount of a fluorescent dye used, the experiment was car- 
ried out in a micro-fluidic set-up based on a recently de- 
veloped technology of soft lithography [14,15]. Moreover, 
we demonstrated just recently that the elastic turbulence 
and effective mixing can bo obtained in a micro-channel 
[16] similar to those in a macro-channel [7]. 

In this Letter we show that a micro-channel is the most 
appropriate system to quantitatively test the recent the- 
ory of the Batchelor decay mixing regime in a bounded 
system. As the result we were able to find the dependence 
of the mixing length on Pe, the critical characteristic of 
a mixer. 

The experiment was carried out in a curvilinear chan- 
nel (inset in Fig.l) scaled down about 30 times, com- 
pared to our macro-channel version [7]. The channel was 
fabricated using a soft lithography [15] out of two com- 
ponents silicon elastomer, PDMS [16]. As a working fluid 
we used a solution of 80 ppm by weight high molecular 
weight polyacrylamide (M^=18 Mda, supplied by Poly- 
sciences). As a Newtonian solvent a 64.4% sugar solu- 
tion in a water with 1% NaCl was used. The viscosity 
of the solvent rjs = 113.8 mPas and the polymer solution 
T] = 137.6 mPas was measured at 22 °C. The polymer 
relaxation time was A ~ 1.08 sec, and the solution den- 
sity was p = 1.31 g/1. The polymer solution was injected 
through two inlets of the micro-channel, and 280 ppm of 
a fluoresccin-labellcd dextran was added in one of them. 
The injection pressure was produced and regulated hy- 
drostatically using vertical rails with rulers and sliding 
stages. We have also checked in successive mixing exper- 
iments performed with the same solutions that the on- 
set of the elastic instability remained unaltered, and the 
flow velocities were reproducible. This suggests that nei- 
ther the viscosity nor the characteristic relaxation time 
changed during experiments. 

One of the critical requirements of the experiment was 
to scan Pe in as wide as possible range of values. We 
varied Pe about 40 times by using tracers with different 
D and also by varying the flow velocity. The former 
has been achieved by using various molecular weights 
fluorescein-conjugated dextran (FITCD) that does not 
contribute to the elasticity of the solution, since it is 
highly branched and stiff. By varying the molecular 
weight of FITCD from lOKda till 2Mda we were able 



to change D of the tracer more than 12 times, from 
5.4 • 10-1° tin 6.6 • 10-9 cmVsec. The values of D in 
water were taken from Ref. [17] and then corrected ac- 
cording to viscosity of a solvent used [11]. 

In order to find a transition to the elastic turbulence in 
a micro-channel both components of the velocity (longi- 
tudinal and radial)were measured by microscopic particle 
image velocimetry (micro-PIV) technique. We used the 
advantage of a narrow depth of a focus of a high mag- 
nification objective (4/im) and an additional pre-filtering 
of the defocused image features. We imaged 0.2 /xm flu- 
orescent particles with a CCD camera mounted on an 
epi-fluorcscent inverted microscope and digitized with 12 
bit, 640 X 512 pixels resolution. As a result, the velocity 
field was measured for each pressure drop with a spatial 
horizontal resolution of 20 x 20 /im and by using a stan- 
dard cross-correlation algorithm between two consecu- 
tive snapshots. An analog of the flow resistance, Aps/V, 
versus the pressure drop per segment, Aps is shown in 
Fig.l. Transition from a laminar regime to the elastic 
turbulence occurs at Aps ~ 50Pa. The average radial 
gradient of the longitudinal velocity, (dV/dr) estimated 
from micro-PIV measurements, allows us to get the Weis- 
senberg number at the transition, Wic ~ 2.7. 

We studied the mixing dynamics in the elastic tur- 
bulence regime along the micro-channel by using a con- 
focal setup (Fluoview FV500 by Olympus) through an 
40 X infinity corrected objective with a numerical aper- 
ture N. A. =0.85. Fluorescent light intensity profiles along 
a single line across the channel as a function of time were 
recorded by a 12 bit photomultiplier with a resolution of 

0. 18 /im/pixel at the acquisition rate of 56 lines/s. 

The interface separating the fluids with and without 
the dye was smooth smooth and parallel to the direction 
of the flow in the case of a solvent injected, although 
little smeared by diffusion for A'' > 45. Whereas in the 
case of a polymer solution injected, the interface became 
unstable at Aps « 50 Pa. First, the interface became 
unstable towards the end of the channel, and then in the 
range of the flow rates used in the experiment on mixing, 

1. e. about twice of the onset value, the flow in the entire 
channel became unstable and chaotic. 

We investigated the evolution of the mixing patterns 
along the channel at various constant flow rates above 
the transition to the elastic turbulence. The maximal 
longitudinal velocity was V ~ 456/im • /sec, and the 
radial velocity component was about 100/im [16]. The 
corresponding average shear rate was obtained from the 
average longitudinal velocity gradient, j'av — {dV/dr) = 
3.7 sec~^ that lead to Re = pVavdi]^^ — 1.6 • 10-^ and 
Wi « 4.0. 

The smallness of Re lead to a conclusion that the mix- 
ing observed was solely due to the elastic turbulence sim- 
ilar to what occurred in a macro-channel [7,10], and so 
we expect to flnd the Batchelor regime of mixing like in 
a macro-channel [7]. We varied Pe more than 37 times 
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from 1.03 x lO" up to 4.14 x 10^. As a result the Batch- 
elor scale, ijb = d - Pe~^^^, changed from 1 x 10"'* down 
to 0.15 X 10~^ cm that was of the same order as a spa- 
tial resolution in the confocal measurements (better than 
0.3/xm). 

As in a macro-channel the mixing in the micro-channel 

in the clastic turbulence regime turned to be rather effi- 
cient as demonstrated in the space-time plots in Ref. [16]. 
A quantitative measure of the mixing efficiency is given 
by the i*'' order moments, Mi, of a probability distribu- 
tion function (PDF) of FITCD concentration, c, defined 
by Mi = (|c— c|') /c*, where c is the average concentration 
of FITCD. We present in Fig. 2 the dependence of Mi and 
M2 on the position along the channel, N, for FITCD with 
molecular weight 70 Kda (at Pe = 5.4 x 10^). In the inset 
of Fig. 2 we show also the decay exponents of the higher 
order moments, 7i, as a function of their number, i, for 
five different Pe values. One notices that the saturation 
occurs for all curves at about the same value of i = 6 
but the saturated value of the decay exponents depends 
on Pe. One can also notice from Fig.2 that, first, both 
moments exhibit exponential decay up to A'' « 40 as pre- 
dicted by theory [5,9]. Second, the decay rate of M2 is 
twice higher as of Mi. And third. Mi is reduced more 
than 30 times and reaches the noise level at A'' > 40, as 
well as M2. Such measurements have been conducted for 
six different values of Pe. The mixing length, measured 
in the number of units, Nmix, and obtained from the ex- 
ponential decay of M2, is presented in Fig. 3 as a function 
of Pe. The data are fitted by L^ix oc Pe^-^^^O-Oi ^j^g^^ jg 
in a good agreement with the theory [9]. 

The latter result indicates that the efficiency of the tur- 
bulent mixing is reduced due to finite size effects, which 
increase the mixing length in an bounded system com- 
pared with that in an unbounded one. In order to estab- 
lish the source of the algebraic dependence of Lmix on 
Pe instead of ln{Pe), as predicted for an imboimded sys- 
tem, we studied the Pe dependence of the mixing layer 
width, W. The latter was identified from the spatial 
distribution across the channel of the first and second 
moments of PDF of FITCD concentration fluctuations 
(see Fig. 4). We associated well-defined peaks in a spa- 
tial dependence of the moments with the boundary layer 
edge (note asymmetry in left and right peak locations 
due to flow asymmetry). The average value of the right 
peak location for the first moment, averaged along the 
channel as a function of Pe is presented in the inset of 
Fig.4 together with the fit: W oc Pe-0-28±o.06_ -jj^jg 

suit is once more in a fair agreement with the theoretical 
prediction W oc Pe~^^^ [9]. 

And finally, the spatial correlation across the channel 
and temporal correlation functions of the FITCD concen- 
tration fluctuations in the bulk and in the boundary layer 
were studied at different Pe. We found that the spatial 
correlation functions decay logarithmically at distances 
above the diffusion length (that is varied from 1 to 4 



/im) and are almost independent of Pe (Fig. 5). The rel- 
atively narrow spatial range of the correlation function 
decay does not allow to reliably distinguish between the 
logarithmic and the power-law behavior, if the exponent 
is small, as predicted by the theory [9] ( the exponent 
~ Pe~^/^). This closeness to the logarithmic decay can 
be explained by the fact that the mixing in the bulk oc- 
curs rather efficiently. The mixing boundary layer per- 
manently supplies the passive scalar into the bulk that 
mimics the Batchclor stationary regime of mixing [8] and 
leads to the logarithmic decay [18]. The temporal cor- 
relation functions in the bulk and in the boundary layer 
(plots are not shown) show similar behavior with about 
twice larger correlation time in the boundary layer than 
in the bulk {ti/2 ~ 13 sec instead of 6 sec in the bulk). 
On the other hand, using the Taylor hypothesis we found 
that the correlation lengths along the the flow are rather 
close in the bulk and the boundary layer {di/2 ~ 2.7 mm 
vs 2.0 mm) and about 200 times larger than the correla- 
tion length across the flow. 
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FIG. 2. The first (squares) and second (dots) order mo- 
ments decay vs. position, N, at Pe = 5.4 x lO''; solid lines are 
the fit. The inset: The decay exponents of the moments, 7i, 
as a function of the moment rmmbor, i, at difl^crcnt Pe: dia- 
monds- Pe = 1.03 X 10*, up-trianglcs- Pe — 2.0 x 10*, squares- 
Pe = 5.4 x 10*, circles- Pe = 1.3 x 10^, down-triangles- 
Pe = 2.5 X 10®. 

FIG. 3. The mixing length, measured in the number of 
turns, Njnix, vs. Pe: squares are the data, solid line is the fit, 
and dotted line shows the logarithmic dependence. 



FIG. 1. Flow resistance, Aps/V, vs. Aps. Arrow indicates 
the onset of the elastic instability. The inset: Schematic dia- 
gram of a curvihncar micro-channel consisting of 64 base el- 
ements smoothly interconnected. The number N of the unit 
counted from the inlet is a natural linear coordinate along the 
channel. 



FIG. 4. Spatial dependence of the first (gray) and second 
(black) order moments of FITCD concentration fluctuations 
across the micro-channel. The inset: Width of the boundary 
layer, W, vs. Pe. Solid line is the fit. 



FIG. 5. Correlation coefficients for the tracer concentration 
vs. Ax/d. Symbols are (1) Pe = 2.5 xlO^ (2) Pe = 1.0 xlO^ 
(3) Pe = 5.4 x 10*; (4) Pe = 2.0 x 10*; (5) Pe = 1.03 x 10*. 
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